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ABSTRACT 


This is a continuation of the work of assembling and testing 
the equipment necessary to conduct gamma-gamma angular correlation 
experiments, started by Fort A. Verser, Jr. at the JUS Naval Post- 
graduate School (1). 

No substantial changes were made in the equipmenr, but rather 
a rerun of the Cobalt-60 experiment was made to confirm the results 
obtained in the previous work. A comparison of the previous work 
confirms that there is a definite distortion in the shape of the 
observed correlation function. 

The statistical analysis of the data and the solid angle 
corrections required for evaluation of the results have been pro- 
grammed for calculation by use of the 1604 Digital! Competer. This 
paper outlines these programs in detail. 

The author wishes to express his appreciation for the assiste 
ance given him by Professor Harry E. Handler and for his guidance 
throughout this investigation. Appreciation is expressed to 
Professor Edmund A. Milne for his advice and comments. Captain T. 
R. Abernathy, USMC, was very helpful in providing valuable assist- 
ance and discussions on various aspects of the computer programing. 
The cooperation and assistance of Mr. Mervyn C. Brillhart in main- 


taining and testing the electronic equipment was indispensable. 
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CHAPTER I 
INTRODUCTION 


The decay of Cobalt-60 to an excited state of Nicke!+60 and 
the angular correlation of the resulting gamma-gamma cascade have 
been studied extensively experimentally (2,6,8,9, and others) and 
the results of the correlation measurements agree very well with 
the theoretical correlation for the decay scheme shown in Figure l. 
Consequently, the decay of Cobalt=-60 is now used as one standard 
for the testing of correlation equipment. 

Angular correlation equipment was assembled at the US Naval 
Postgraduate School by Verser in 1960 for the investigation of the 
Properties of low-lying excited nuclear states. His testing of the 
operating characteristics of the equipment by the angular corre- 
lation of the cascading gamma rays resulting from the decay of 
Cobalt~60 yielded an observed correlation curve whose shape was not 
in agreement with the theoretical curve. This work was initiated 
in order to attempt to verify Verser'’s results. 

In order to facilitate the statistical analysis of the data, 

a program has been worked out for the use of the 1604 Digital 
Computer. The program follows the procedures outlined by M.E. 
Rose (2) and is contained completely in Appendix 1. 

An additional program has been written to compute the solid 
angle corrections, which must be applied to the theoretical curve, 
by the method of M.E. Rose (2). This progam is outlined in Ap- 
pendix IT. 

The results of this work agreed within statistical expectation 
with those obtained by Verser, indicating that the unexpected shape 
of the curve is apparently real and not an instrumental effect. 


These results are given in Chapter VI. 





PARTIAL DECAY SCHEME OF COBALT-60 


COBALT-60. 
5 4 ; 2.814 MEV 
NICKEL-60 | 
4+ —> ~2.505MEV | 
2 + _ : 1.333 MEV 
| O + S O MEV 
SPIN&PARITY . _ ENERGY 


FIGURE | 2. 





CHAPTER II 
THEORETICAL CONSIDERATIONS 


The determination of the angular momentum quantum numbers, 
or spins, of low-lying short-lived excited nuclear states by 
the nuclear spectroscopic tool, known as angular correlation of 
successively emitted radiations, is quite common practice today. 
Only within recent years, because of the vast advancement of 
electronics, has practical utilization of this tool become 
possible. 

Since nuclei, under ordinary conditions, are randomly ori- 
ented in space, it is not possible to observe a radiation pattern. 
This is due to the fact that the probability of emission of a 
radiation by an excited nucleus, depends on the angle between the 
nuclear spin axis and the direction of the emission. Therefore, 
in order to observe a radiation pattern, it is first necessary to 
orient the nuclei. 

One method of orienting certain nuclei consists of placing the 
sample in an electric field gradient or a strong magnetic field at 
a very low temperature. Another method, the one used in the present 
consideration, consists of selecting nuclei whose nuclear states 
have spins lying in preferred directions. This happens to be the 
case for the intermediate state if a nucleus decays bv successive 
emission of two radiations. 

In an angular correlation experiment the first radiation is 
observed in a fixed direction. This establishes an axis to which 
the direction of emission of the second radiation, originating from 
the state formed by the first, can be referred. The second radiation 
has preferred angles of emission with respect to the direction of 
the first. 

In order that the correlation exhibit maximum anisotropy, the 


angular momentum vector of the intermediate state must not change 





direction significantly before the second radiation occurs. The mean 
life of the intermediate state must be less than about ae seconds, 
the typical precession period of the nuclear spin about local per- 
Eurbing fields. The angular correlation function is determined by 
the nature of the radiations and by the spins of the states involved, 
and thus its measurement may lead to an unambiguous choice of spins 
for the states. 

The theoretical expressions for the correlation functions have 
been worked out for a large number of cases of interest (7). For 
the particular case of a gamma-gamma cascade, the correlation 


function W(@) is 


n n 
max max 

WO) = > eee a(Cosmaye = > B Cos'0 (1) 
n = Q n = 0 


where PY is the Legendre Polynomial of even order n, Dax is the 
smallest even integer of the set of three numbers consisting of 
twice the multipole order of each of the two gamma rays, and 9 is 
the angle between the directions of emission of the two gamma rays. 
The theoretical correlation function for the cascade in Nickel-60, 


assuming point detectors and a point source, is 
W(8) = 1 + 0.1020 P,, (Cos Q@) + 0.0091 P, (Cos 0) (2) 


with an anisotropy 


g = W(180") - W(90") = 91667 (3) 
w(90°) 


Many workers, after correcting these theoretical results for the 
effect of non-point detectors have found good agreement between them 


and their experimental correlation results. 














CHAPTER III 
DESCRIPTION OF EQUIPMENT 


The equipment needed for the gamma-gamma correlation 
measurements must furnish means for detecting the radiation, 
pulse-height amplification with selection or discrimination, 
coincidence analysis, and adequate recording of the counts in 
the various counting channels. A mechanical system must also be 
provided for mounting the detectors and changing their relative 
positions and for holding and positioning the radioactive source. 

The apparatus used for this work was that assembled by Verser 
and described in detail in his report. Therefore, only a brief 
description of the equipment and of the minor alterations per- 
formed on it will be presented. 

Figure 2 is a photograph of the mechanical setup. The origi- 
nal base of the correlation table has been replaced, the new base 
is of solid brass with a bearing mounted shaft to enable accurate 
manual detector positioning. 

The electronic equipment (Figure 3) has the functions of 
detection of the radiation, pulse-height amplification and dis= 
crimination, coincidence analysis, and recording the number of 
counts. The only change in the electronic equipment was in the 
pulse-shaping networks used in the initial calibration of the 
apparatus. The pulse shapers were rebuilt to give a better simu- 
lation between the test pulses and the actual pulses encountered 


under operating conditions. 
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CHAPTER IV 
EXPERIMENTAL ANALYSIS 


The angular correlation function W(@) (Equation 1) must be 
expressed in terms of experimentally-observable quantities. The 
equipment records the total number of coincidences as a function of 
the angle between the detectors, and the total number of radiations 
captured in each channel. Ch. the total number of coincidences in 


a counting period, is given by 


Ge = CO. Wee Ce eG 4 
t g a b SoD, 
where 
C; = the mimber sor cenugine coincrdencessduerto thevtwe 
6 cascading radiations, 
oe = the number of accidental coincidences due to the 
finite resolving time of the electronic equipment, 
and 
C, = the number of background coincidences. 


The half-life of the source is long compared to the time required 
to collect the necessary data; therefore, equation (4) can be re- 
arranged and rewritten in terms of counting rates. 


C= 1G = o =aG (5) 


g t b 


S. can also be expressed in terms of other parameters: 


oe oe ‘An W ( 
Co CE Pay MaaP on + Epa 2% 10% 1% 2180179428 ag OP i) (6) 
where 
N = the absolute source strength 
P ik = the fraction of nuclei decaying by the emission of 


gamma ray "j'’ and gamma ray "k"' in a cascade 





ij 


ij 


Ores 
LJ 


WO) 


the fractional solid angle subtended by the crystal 


in Channel “a1 


the fraction of gamma rays "'j'' which leave the source 
in the right direction to enter crystal "i" that are 
not absorbed before reaching the crystal 

the total verfhicrency of the crystal in Channel "i 


for gamma ray "j" 


the fraction of the number of gamma rays "j"’ captured 
by the crystal in Channel "i'' whose pulse heights are 
accepted by the energy discrimination of Channel "i" 
the angular correlation function at an angle ee 


averaged over the solid angles used. 


N., the counting rate in Channel "i'', can also be expressed as a 
i 


function of certain of these same parameters. 


where 


ib 


i 


il 


i 


a, 95%) 
i Pi NoM a eo (7) 


the number of gamma rays "i'' per disintegration 


the background counting rate in channel "i" 


From Equations (5) and (6) we get for the correlation function, 





(8) 


where f is the appropriate function of the €'s, Q's, &'s, and G's. 
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By combining Equations (7) and (8), we get for the correlation 


function 


(C,-C,-C.) Ne 





(9) 


wOevere vised, function of there cy Gs. © seis 60 o. 


In this form we) is relatively insensitive to equipment stability. 
Therefore, it is only necessary to calculate the modified correlation 


function given by 





= ey aay 
w'(@) = - 


x ee =aaeae (10) 
7 (NWN? Na N2b? | H 








The correlation runs were done with the movable detector at 
seven angles: 90°, 105°, 120°, 135°, 150°, 165, and 180°. A 
minimum of six two~hour counting periods were completed at each 
angle, and the angular order in which they were taken was varied 
to minimize any systematic fluctuations. Each run yielded single 
channel counting rates of approximately 1200 counts per second and 
a total coincidence rate of approximately 1.1 counts per second. 

Before the experimental results can be compared to theory, the 
data must be checked and corrected, if necessary, for the various 
deviations from the ideal arrangement, i.e., centered point source, 
no spurious coincidences due to the coincidence resolving time, no 
Scattering, no background, and perfect stability of the electronic 
equipment. In addition, the theoretical correlation function must 
be corrected for non-zero solid angle detectors. The methods for 


accounting for these effects are discussed in the following para- 


graphs. 


il 





A. SOURCE SIZE AND POSITIONING 


As a test of the positioning of the source in the center of the 
correlation table, the integral counting rate of a portion of the 
Cobalt-60 spectrum was measured in each channel as a function of the 
angle between the counters. This was done by Verser and checked again 
after the base of the correlation table had been changed. Within sta- 
tistical expectations, the counting rates were constant. The source has 
a volume of about 0.001 cubic inches and can be regarded as a4 poink source 


at the distance at which the detectors were placed from it. 


B. SCATTERING 


Compton scattering occurring outside the detectors can give rise to 
unwanted coincidences. These in general will tend to smear out the mea- 
sured correlation function. This problem is minimized by use of scintil- 
lation crystals as detectors and accepting only the full-energy peaks of 
the desired gamma rays. The method for determining the pulse height se- 


lection limits has been discussed in detail by Verser. 


C. BACKGROUND CORRECTION 


With the discriminators at the same settings as were used during the 
correlation runs, background runs were made before and after the corre- 
lation runs. The background coincidence counting rate was negligibly small. 
and the background rate of each channel was such as to yield total backe 
ground counts of the order of 0.4 percent of the total counts in each 
channel during the correlation runs. These background rates were used as 


corrections to the gross counting rates obtained during correlation runs. 


D. DETERMINATION OF GENUINE COINCIDENCE RATE 


Since the background coincidence counting rate proved experimentally 


12 





to be negligible, Equation (5) becomes 


° 


For the total counting time, T, the total number of coincidences is Ct 


and the standard deviation in the total coincidence rate is 


ae || ees Ci) 


Therefore, the standard deviation of the genuine coincidence rate is 


C C 
a 


a ee oat (12 
g t a i T 





The accidental coincidence rate is given by 


o 9 9 ° D 

©. i am a ot: e (13) 
where 

TY = Coincidence resolving time of the coincidence analyzer. 


in order to minimize the accidental coincidence rate, the source 
strength must be kept as small as possible. However, if the source 
strength is decreased, the counting time must be increased to obtain a 
specified precision in the counting rate. Thus, a compromise must be 
reached between precision and counting time. 

The smallest resolving time obtainable with the present equipment is 


experimentally determined to be 
Ve = (0.1820 + 0.0016) microseconds; 


this resulted in an accidental coincidence rate that was about 390% of the 


13 





total coincidence rate. The approximate source strength of the water 
solution of CoCl, used in this experiment was 0.05 millicuries of Cobalt-69, 


E. REFERENCE COUNTING RATE 


With the discriminator in each channel adjusted to accept only the 
full-energy peaks of the desired gamma rays, the integra! counting rate 
of the Cobalt-60 spectrum in each channel was determined. This gave a 
reference counting rate for each channel. During the course of the corre- 
lation runs, either the discriminator settings or the gains of the linear 
amplifiers were adjusted before each run so that, the integral counting 


rate agreed within + 5 cps with these reference rates (Table 1). 


TABLE 1 
REFERENCE COUNTING RATES 
(1065. 2654 0257) cps 





A TS Ce Bs tre re 


Fas e 
I 


Ne pC 7322 0.59) weeps 
eee AS eee a is: ee A a en 2s es biel) gS 
As a measure of the gain fluctuations during each counting run, the 


integral counting rate of each channel was measured after the run using 
the same discriminator setting that was employed during the run. The 
difference between this counting rate and the reference counting for that 
channel was determined. If the variation exceeded + 20 cps in either 
channel, the run was discarded. 

This procedure resulted in about 30% of the runs being discarded 
because they fell outside of the maximum acceptable fluctuation. This 
indicates that the electronic system is not as stable as it should be. A 
check of the line voltage revealed substantial fluctuations during any 24- 
hour period. The temperature of the room was not controlled in any manner 
and varied considerably. These could have been the major causes of the 
large counting rate fluctuations since the discriminators are highly sensi- 


tive to both temperature and voltage variations. 
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F. SOLID ANGLE CORRECTIONS 


The solid angle correction which must be applied to the theoretical 
correlation function to account for the non-zero solid angle of the de- 
tectors and enable a comparison to be made with the measured function is 
described by Rose (2). The correction for two cylindrical crystals whose 
axes intersect at the source involves the numerical evaluation of integrals 


of the form 


4 TX.) 
lh akon P (Cos @) ;1 -e ‘| Sin @d@ (14) 
cee | n C 2 
=o 
where 
ie = the full-energy absorption coefficient of the detector in 
Channel "i" for the gamma ray; 
o = the half-angle subtended by the front face of the crystal 
ieowannel. a; 
= ai 
X.(@> = t.Sec@ for 0<G= Tan : ee 
: = h. +t. 


-1 r 
= =- < 
x. (9) r .Cscf h,Sec@ for Tan hoa? f : 


eee 
i it 
where 
n. = the distance from the source to the crystal in Channel "i"'; 
t. = the thickness of the crystal in Channel "i"'" 
r. = the radius of the crystal in Channel "i". 


The attenuation factors, Q. are functions of these integrals: 





Q =a : tes) Mz (lap 
Oxi One 





—_— 


The corrected correlation function becomes 


7 Nmax 
W(0) = Q An en 6°89) 
n=o 


i 





CHAPTER V 
EXPERIMENTAL RESULTS 


The modified correlation function for the nen run at angle O- given 
by Equation (10) is 
7 (0.) = oe = Ce = C. 
ee mn i a 
(N7-Nj,) O-N>,? 
k 


This means that the seven counting rates es C.> Ch Ni» ND? Ny» and Nop 
must be determined for each run at each angle and the standard deviation 
of each rate calculated. 

The background counting rates Ni and Noy for each run were experi- 
mentally determined. The single-channel counting rates corrected for 
scaling losses yielded N, and N, The accidental counting rates were cal- 
culated for each run from Equation (13). The total coincidence counting 
rate, ae was experimentally determined. The background coincidence rate, 


Cy» was negligible throughout the experiment. 

The experimental value of the modified correlation function We (9. ) 
and its standard deviation were calculated for each run at each angle. 
From these the weighted mean W' (9. ) and its standard deviation were com=- 
puted for each angle. 

The weighted mean of the measurements taken at each angie was deter- 
mined by using the inverse square of the standard deviation of each mea= 


surement as a weighting factor: 


Ze 
6} (18) 








where 
Ne = number of runs at the angle O.. 
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The standard deviation of the weighted mean (W') was calculated from 


the individual standard deviations: 


ee / (19) 


A comparison between the standard deviations of the means calculated 


from Equation (19) and those calculated from 


lw’) = | eae (20) 


revealed no significant differences. Therefore the standard deviation 
obtained from Equation (19) were used in subsequent calculations and the 
method discussed in paragraph E Chapter IV for treating gain instabilities 
is realistic. 

The values of the weighted means and their standard deviations are 


Rapulacted =inelable 2. 


TABLE 2 
The Weighted Means w' (0. ) 





0, (degrees) Weighted Standard 


Me an Deviation 











ae 





oe 


l= 


mental wie ) by the method of least squares as outlined by Rose (2). 


A curve of the form w'(@) = AL P (Cos®) was fitted to the experi 
This analysis was programmed to be carried out by the 1604 computer. (See 
Appendix II for details). 


The values of At are given by the matrix equation 


fed (21) 








H aT al at lat 





e 


a 


| a’ 


where 
| Al| = the 7 x 3 matrix consisting of the Legendre Polynomials, 


Po Po» and P,, evaluated at the seven angles (Figure 5A). 


4? 


tH 


|] a7 


| wll - a diagonal matrix with elements consisting of the inverse 


the transposed A matrix 


squares of the standard deviations of the experimental 


w' (0. ) (Figure 5B) 


a one column matrix whose elements are the experimencal 


values of w'(0.) (Figure 5C) 


| £2 


i 


| al = the one column matrix whose elements are respectively 


the least-square coefficients A‘ A, and A; (Figure 5D) 


=] : : : , . 
|| aTuall = the inverse of the matrix obtained by multiplying the 
three designated matrices together. The diagonal terms 
of this matrix are the squares of the standard deviations 


of the least-square coefficients, A. 


The experimental values for the coefficients A’ and their standard 
deviations, normalized so that A’ = 1, are listed in Table 3. The corre- 


sponding values from Verser's work are included for comparison. 
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TABLE 3 


Least Square Coefficients 


Verser This Work 





L9 





MATRIX CONF IGURATIONS 


| 
P (Cos90”) P,(Cos90°) P,, (Cos90°) 
P_ (Cosl05") P,(Cos105") P,(Cos105°) 
: O O fee) 
P ,(Cosi20 ) P, (Cos120 ) P, (Cos120 ) 
| P (Cos135°) P, (Cos135°) P, (Cos135°) 
| O O O 
| P (Cos 150 ) P.(Cos150 ) P,, (Cos 150 ) 
| P (Cos165 ) P,(Cos165. ) P,(Cosl65-) 
; P_(Cos180~) P,, (Cos180°) P, (Cos180~) | 
(A) 
O-*(90°) 0 0 0 0 0 0 
0 Oa-cles no 0 0 0 0 
0 0 O~* (120°) 0 0 0 0 
0 0 0 6245 0 0 0 
0 0 0 0 Ca (isds 6 0 
0 0 0 0 0 O65 0 
0 0 0 0 0 0 T2180) 
(B) 
w' (90°) 
W' (105°) 
W' (120°) A 
(135 4) A 
w'(150°) A, 
W' (165°) 
w' (180°) ep) 
(C) 
FIGURE 5 
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A. Experimental Anisotropies 


The anisotropy is defined: 


R —- Ww(180°) - W(90°) | (22) 
w(90°) 


The anisotropy for the least~square curve was calculated from the experi- 
mental correlation function. Since the anisotropy can be cdtermined with 
much greater precision than can the least-square coefficients themselves, 
1t provides a more sensitive means for comparing the experimental and theo- 


retical functions. The experimental anisotropy is 
R = 0.167 + 0.013 


The standard deviation was calculated by the formula presented by Klema 
and McGowan (6). The computer program for the calculation of R and its 


standard deviation is outlined in Appendix I. 





CHAPTER VI 
RESULTS AND CONCLUSIONS 


A. COMPARISON BETWEEN EXPERIMENT AND THEORY 


The experimental and corrected theoretical functions for this experi- 


ment are: 
Experimental: 


W'(@) = (1.0000 + 0.0026) + (0.0971 + 0.0057) P, 
+ (0.0231 + 0.0066) P, 


Theoretical: 


W(9) = 1.0000 + 0.0996 P, + 0.00838 Be 


2 
The experimental anisotropy was calculated to be Re Ores + 0.013 and 
is compared with Verser's value on es Om + 0.027 and with the theo- 


retical anisotropy corrected for the solid angle, R = 0 oie 


Figure 6 shows the experimental results and their standard deviations for 
this work and that of Verser. Also shown are the least=square correlation 
curves for both sets of data and the corrected theoretical correlation 


ae : O 
curve. All the curves have been normalized to unity at 9 = 90 , 


B. DISCUSSLON AND CONCLUSIONS 


From Figure 6, it can be seen that the curve obtained in this work 
agrees with that obtained by Verser. Both of these curves show an ap- 
parently real dip in the region from 105° to 135°. It is only a remote 
possibility that this dip can be attributed to statistical fluctuations, 
and it is unlikely to be due to instrumental effects. 

Two conjectures of possible mechanisms are; 

(a) The source may have crystallized resulting in preferential 
alignment of the nuclei within the crystal due to internal 


electric, fields. 


a2 
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(b) The ferromagnetic properties of Cobalt could have led to 
a preferential alignment of the nuclei in the earth's mag- 
HEtrelfield:. 

The following series of experiments are suggested for future work. 

(a) The results for measurements conducted for angles between 
180° and 270° should be compared with those already ob- 
tained. 

(b) In order to check on the possible effects of orientation 
due to crystalline electric fields, the experiment should 
be repeated with the present source at various orientations 
relative to the fixed counter. 

(c) In order to check on the possibility of orientation by mag- 
netic effects the source should be placed in a magnetic field 
of known direction and the measurements repeated. 

(d) A correlation curve using a powdered Cobalt-60 source in 
place of the present source should be obtained for comparison. 

These experiments would determine conclusively if any effects causing pre~- 
ferential nuclear orientation are present. 

The following instrumentation improvements are suggested in order to 
improve the stability of the equipment. 

(a) A very stable power supply should be purchased to aid in 
preventing the line voltage fluctuations from interfering 
with the electronic system. 

(b) The effect of temperature variations on the electronic system 
can be minimized by placing the equipment in a temperature 
controlled room for conducting future work. 

(c) A slow coincidence unit whose resolving time is independent. 
of the counting rate should be obtained. 

The precision of the results should be improved by the addition of a 
fast coincidence channel working in conjunction with the slow coincidence 


unit. to reduce the number of accidental coincidences. 
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Since the experimental curves using Cobalt-60 have not agreed with 
the theoretical curve, the experiments on Cobalt-60 should be continued 
until the source of the distrotion is ascertained. Perhaps the effect 


can be related to the physical properties of the coCl, source, 





APPENDIX I 


Least-Square Fit of Experimental Results 


The comparison of the experimental correlation function to the 
functional form of the theoretical correlation function given by Equation 
(1) is accomplished by the method of least-squares, as outlined by Rose 
(2). The analysis requires the solution of Equation (20) for the values 
of A’: 

The computer solution of the matrix equation, like the solid angle 
calculation, is carried out in floating point format. In addition to the 
solution for the values of Ai. the program also calculates the anisotropy 
by Equation (21), and its standard deviation by the formula presented by 
Klema and McGowan (6). A flow sheet outline of the program is given in 
Figure (7). In order to provide a certain amount of flexibility to the 
program, the angles used have been left as a parameter to be supplied. 
The program as written is, however, limited to seven different angular 
positions. The parameters needed for the program may be placed anywhere 
in the computer and their addresses supplied to the program by placing 
them in the B registers in the following order: 

Bl = weighted means of the correlation function in ascending angular 

order. 

B2 - the standard deveations of the weighted means in the same order 

as above. 


B3 = the values of the angles used in ascending order. 


The program is available for future use on both bioctal tape in 
machine language, and IBM cards in assembly language. The bioctal tape 
contains the following: 


1. The main program for calculation of A’ 
2. The program of normalization of the AL 


The program for calculation of the resulting correlation function 
4. The normalization of the correlation function 


5. The calculation of the anisotropy 
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6. The calculation of the standard deviation for the anisotropy 


i= the» subroutines 


a. 


C. 


Floating Point Sine-Cosine 
b. Floating Point Square Root 


Floating Point Decimal Output. 


The results of the computations are dumped on magnetic tape unit (4) 


in the following sequence: 


Lines | thru 3 


Lines’ 4 thru 6 


Line 7 


Line 8 


Line -9 thru 13 


Line 


Line 


Line 


Line 


Line 


14 


16 


ry 


18 
ee, 


and 


L5 


The matrix inverse 


Check answer; included in the program is a check 


of the inverse which consists of multiplying the 


inverse by the original matrix 


The values of Ao» A anaes A 


The 


The 
The 
The 


The 


The 
The 


Z 4 


standard deviations of the A’ 


resulting correlation functions 
noLtMmalrzed correlation functions 


normalized An 
normalized standard deviations of the An 


anisotropy 


square of the standard deviation of the anisotropy 
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STATISTICAL ANALYSIS PROGRAM 


ORG 12000 
REM BURTON J CONWAY APRIL 1961 
REM LEAST SQUARE FIT OF ANGULAR 
REM CORRELATION DATA. 
REM §$ MUST BE SUPPLIED WITH 
REM §§ PARAMETER ADDRESSES AS FOLLOWS. 
REM Bl CORRELATION FUNCTION. 
REM  B2 SIGMA OF CORRELATION FUNCT. ON. 
REM B3 ANGLES USED IN DEGREES. 

ZEUS EQU 13000 

RHO EQU = = 13100 

TRANS EQU 13200 

ATW EQU 13230 

ATWA EQU 13260 

INVER EQ 13500 

CHECK EQU 13520 

IAT EQU 13540 

TATW EQU 13570 

ANS EQU 13620 

FINIS EQU 13630 

SINFL EQ 60000 


MATRIX EQU 1.0000 


a 





12000 


12001 


12002 


12003 


12004 


12005 


12006 


12007 


12010 


12011 


L201e 


56 oL 
50 1 
at 2 
56 2 
50 2 
50 0 
56 3 
Do, 
ee 
32 0 


15 4 


50 0 
50 0 
50 O 
50 0 
201 
50 0 
54 1 
1? 0 
121 


pow 


12045 
00000 
12032 
12033 
00000 
00000 
12004 
00000 
00000 
12512 
60001 
00000 
00000 
00000 
00000 
00000 
124.63 
00000 
00006 
12004 
12463 
12463 


ANISO 


START 


NEXT 


ECU 
SIU 
ENT 
onl 
SIU 
ENT 
ENT 
SIU 
ENT 
LDA 
FMU 
SLJ 
ENT 
ENI 


ENT 


ENT 
oLrA 
ENT 
ISK 


SLJ 


30 


PW Ww 


© 


13624 
FORMMU 
O 
SIGMA 


SIGMA+1L 


CONV 
SIN®L+L 
O 


O 


ON 


START 
TieMP STARTS FORMATION 


TEMP OF THE LEG. POL: 





22019 


L2014 


2015 


12016 


T2017 


12020 


i202). 


12022 


12025 


12024 


La0ea> 


12026 


oe 
a 
DD 
20 
L2 


22 


ve 
32 
32 
31 
30 
DD 


20 


20 
D1 
2 


20 


aL 


0 


12h61 
12320 
12316 
12317 
12326 
12461 
12322 
12462 
12461 
12461 
12323 
12462 
12320 
12321 
W220? 
00000 
00006 
12012 
12325 
12434 
OO0O]L 
1.2326 
12434 


00001 


FORM 


FMU 
SB 
KDV 
STA 


LDA 


FSB 
FAD 
FDV 
STA 
ENI 
ISK 
SLI 


LDA 


ora 


INI 


ol 


NEXT 


WORKL 


ABLE, 


WORK2 


ABLE 





me0e [ 


E2050 


12051 


12032 


12053 


12034 


12055 


12036 


22057 


12040 


12041 


12 
20 
Bi 
50 
54 
2 
2 
D2? 
22 
20 
54 
ie 


50 
12 
22 
12 
20 
51 
50 
54 
12 


2 


0 


2 


0 


W252) 
12434 
00001 
00000 
00006 
12024 
Hes 16 
00000 
00000 
12463 
00006 
12032 
00000 
00000 
12555 
12043 
12h63 
1255) 
00001 
00000 
00006 


12036 


SIGMA 


FORMCAD 


STA 
INI 
ENI 
ISK 
oLJ 
LDA 
FDV 


FDV 


ISK 
SLJ 


ENT 


LDA 


LDA 
STA 
INT 
ENI 
ISK 


SLJ 


32 


WORK4Y 


ABLE 


FORM 


FLT 


TEMP 


SIGMA 


FORMCAD 





12042 


12043 


12044 


12045 


12046 


12047 


12050 


12051. 


12052 


1205) 


12054 


L2055 


12 0 


wr 
O 
O 


120h4 
00000 
00001 
12036 
00000 
00000 
00000 
234k 
00006 
12045 
00000 
00000 
L00CO 
00000 
00000 
00000 
00000 
00002 
00007 
12h34 
00003 
06626 
00000 


13200 


FORMMU 


SLJ 


ENT 
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SLJ 
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ENI 
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ENI 
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ENT 


ENT 


ENI 


ENI 


ENT 
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J 
ra) 
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O ABLE 
05 
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= 
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O 
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ENI 


O TRANS 


O CAD 


O ATW 


Ot 
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12077 
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ee 
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50 


O 
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10000 
00000 
00000 
00000 
00000 
00005 
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o0ck4 
13272 
00000 
13500 
10000 
00000 
00000 
00000 
00000 
oo000k 
00003 
13500 
00003 


13260 


ENI 0 3 
ENI O ATWA 
SLJ 4 MATRIX 
ENI OO 

ENI OO 

ENI 00 

ENI 00 

ENI O 5 

ENI 0 3 

ENI O ATWA 
ENI O 36D 
ENI O ATWA+12 
ENI 00 

ENI O INVER 


SLJ 4 MATRIX 


3 
an 
ae 
Oo 
= sS 


2 
bh] 
S 
LN 


ENI O INVER 
ENI O 3 


ENI O ATWA 


Ze 


INVERSE 


CHECK 





121.05 50 O 00003 ENI O 3 


50 0 13520 ENI O CHECK 
12106 75 & 10000 SLJ 4 MATRIX TAT 

50 0 00000 ENI 00 
12107 50 0 0000 ENI OO 

50 0 00000 ENI 00 
12110 50 0 00000 ENI 00 

50 0 00004 ENI o4 
Tan. 50 0 00003 ENI O03 

50 0 13500 ENI O INVER 
12112 50 0 00003 ENI 0 3 

50 0 13200 ENI O TRANS 
12113 50 0 00007 ENI O07 

50 0 13540 ENI O IAT 
12114 75 4 10000 SLI 4& MATRIX TATW 

50 0 00000 ENI 00 
12115 50 O 00000 ENI OO 

50 0 00000 ENI 00 
12116 50 0 00000 ENI OO 

50 0 00004 ENI O 4 
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7 
CAD 


TT 





12134 


12155 


12136 


12157 


12140 


12141 


12142 


12143 


12144 


50 0 
Th 7 
75 4 
50 0 
Ol O 
O4 O 
12 0 
70 0 
54 4 
15 0 
54 1 
tooo 
54 1 
(ee 
50 0 
Th 7 
15 4 
50 0 
Ol O 
ob O 
12 0 


7O O 


00000 
4200 
71000 
00000 
13500 
00001 
12315 
12136 
00002 
12135 
T1777 
12141 
TTTTT 
12142 
00000 
42040 
71000 
00000 
13520 
00004. 
123515 


12145 


ENI 
XE 
SLI 
ENI 
OL 
oly 


LDA 


ISK 


SLJ 


ISK 


SLJ 


ENI 


oLJ 


OL 


LDA 


38 


00 

7 4eoko 
4 71000 
00 

O INVER 
oval 

O MIKE 
0 /-1 
42 

0 /-3 
LL 
0 / 
ce 
0 / 
00 

7 4eoko 
4 71000 
On0 

O CHECK 
o4 

O MIKE 


0 f/-l 





12147 


12150 


12151 


teipe 


12155 


1215) 


12156 


Lolo? 


12161 


sl. h 00002 


75 0 121k4 


S41 T7777 
7500 12150 


50 0 00000 


7h 7 4eoho 


75 4 71000 


50 
OL 
Ou 


re 
od 
De 
30 


O 


O 


Q 


by 
l 


O 


00000 
13620 


00007 


13620 
12434 
13100 
O000L 
13622 
1243h 
13100 
13100 
13622 
00001 
12434 


33100 


Isk 4 2 


SLJ 
ISK 
SLJ 
ENI 
EXF 
SLJ 
ENI 

Ol 

Ou 


REM 


LDA 


INT 


LDA 


FAD 
oT A 
LDA 
INI 
PMU 


FAD 


39 


0 


1 


0 


O 


i= 
TITTT 
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0 


7 42040 
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O 


- oO 
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eae 


(1000 


STARTS FORMATION OF THE 
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ABLE 


RHO 


L 


ANS+L 


ABLE 


RHO 


RHO 


ANS+2 


ABLIS 


RHO 


AOFO IN STORAGE 


SePe TN A REG 
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A REGISTER 


AUPL TN A REG 


SUM DY 4 REG 





12162 


12163 


12164 


12165 


12166 


12167 


nel [0 


peal ad 


ole 


12173 


12174 


20 6 13000 
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Si 4 oocek 


12 
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15 
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12 
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d 
yy 
a 
O 


O 


12154 
13500 
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00000 
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Ok O 00010 
S41 TTTTT 
75 0 12176 
50 0 00000 
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12 1 13000 
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51 1 00003 
54 2 00006 
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50 0 00000 


74 7 “eoko 


ENI 
EXF 


SLI 


Ol 


LDA 
RAD 
ISK 
oLJ 
LDA 
FDV 
STA 
INL 
ISK 
SLJ 
ISK 


oLJ 


XE 


4] 


O TTA 

O 10 

tae 

0 / 

00 

7 4eoko 

4. 71000 

0 0 

QO ZEUS 

O11 

O MIKE 

O /-l 

5 4 

O /-3 

1 ZEUS NORMALIZES THE 
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2 NORM FUNCTION 
see 

2 6 

0 /-2 

TTT 

0 / 

00 

7 42040 





12215 


12214 


12215 


12216 


Le2l] 


12220 


12 
50 
OL 
O4 
12 
70 
54 


19 
50 


4. 71000 
© O0C00 
0 12472 
O 00016 
O 12315 
O 12212 
2 00001 
O 12211 
1 00000 
0 00000 
1 13620 
O 13620 
1 12507 
0 00000 
1 00002 
O 12216 
a ee 
QO Le2el 
OQ 00000 
7 Leoko 
4 71000 


QO 900000 


SLJ 
ENT 
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LDA 


ISK 
SLI 
ENT 
ENT 
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ENT 
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SLJ 
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21 

QO f-3 

10 
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00 

12 
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oy 

0 0 
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122350 
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Tee5c 


Lee5s 
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12259 


12236 


12257 
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ob O 
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D5 0 
20 1 
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13 0 
54 1 
15 0 
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75 4 
50 0 
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ok oO 
12 0 
ico 
33 0 
20 0 
She 
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50 0 
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00020 
12501 
13000 
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00000 
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00000 
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APPENDIX LIL 


Solid Angle Correction 


The solid angle correction, which must be applied to the theoretical 
correlation function to enable a comparison to be made with the experi- 
mental correlation function, is described by Rose (1). The correction 


involves the numerical evaluation of the integrals given by Equation (14): 


1 


ee f P_(Cosd) Cimeeey i) ead a0 
O 

and the subsequent solutions of Equations (15) and (16). 

The computer solution for this correction was carried ont in floating 
point format. The numerical integration was computed using an interval 
of 0.091 radians. A flow sheet outline of the program is given in Figure 
(8). The program must be supplied with the parameters ee ie Eo and st 
same in both channels (see Table 4) thus enabling solid angle corrections 
to be made for any experimental set-up to which these corrections may be 
applicable. These parameters may be placed anywhere in the computer and 
their addresses supplied to the program by placing them in the B registers 
in the following order 

Bl = r (address of the crystal radius in centimeters) 

B2 t (address of the crystal thickness in centimeters) 


B3 


t! 


h (address of the crystal to source distance in centimeters) 
B4 ~ T (the value of the full-energy absorption coefficient of the 


: =] 
detector for the gamma ray in centimeters ") 


The values of the parameters used in the calculation of the corrected 
angular correlation coefficients for the theoretical correlation function 
and the results of the calculations are listed in Table 4. 

The program is available for future use on bioctal tape in machine 
language and on IBM cards in assembly language. The bioctal tape contains 
all the necessary subroutines to enable a complete computation, i.e. 

1. The main program for solution of Equation (14), 


2. The subroutines 
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dk ey Ny rel aa 


a. Floating Point Sine-Cosine 
b. Floating Point Exponential 
c. Floating Point Decimal Output. 
The results of the computations are dumped on magnetic tape unit (4) 


in the following sequence: 


Line 1 Ly 
Line 2 1, 
Line 3 I, 
Line 4 T/T 
Line 5 I,/1o 
Line 6 Q, 
Line 7 AQ, 
Line 8 Q, 
Line 9 A,Q), 
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TABLE 4 


SOLID ANGLE CORRECTION PARAMETERS 


h (cm) 
t (cm) 
r Cem) 


7 (nee) 
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REM BURTON J CONWAY APRIL 1961 
REM SOLID ANGLE CORRECTIONS 

REM OF THE THEORETICAL CORRELATION 
REM CURVE. MUST BE SUPPLIED 

REM §§ WITH PARAMETER ADDRESSES AS 
REM —- FOLLOWS. 

REM Bl ADDRESS OF XTAL RADIUS. 
REM Be ADDRESS OF XTAL THICKNESS. 
REM B83 ADDRESS OF CRYSTAL TO 

REM § SOURCE DISTANCE. 


REM B4 ADDRESS OF PARAMETER TAU. 


HOLD ECU 5500 
TEMP EGU 6500 
WORK2 EQU 7500 


WORK: ECU 10500 


HANS ECU 11500 
SOLD ECU 12500 
JAYO ECU 13500 
JAY2 ECU 13504 
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